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ABSTRACT 


One  of  the  difficulties  in  designing  controllers  for  infinite- dimensional  systems  arises  from 
attempting  to  calculate  a  state  for  the  system.  In  this  paper  it  is  shown  that  Galerkin 
type  approximations  can  be  used  to  design  controllers  which  will  perform  as  designed  when 


implemented  on  the  original  infinite-dimensional  system.  No  assumptions,  other  than  those 
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1.  Introduction 


There  are  computational  difficulties,  apart  from  the  theoretical  problems,  to  designing 
controllers  for  a  class  of  systems  whose  dynamics  are  described  by  partial  differential  equa¬ 
tions.  Unlike  finite-dimensional  linear  dynamical  systems,  which  are  often  easily  solved  for 
the  state  of  the  system,  a  closed  form  solution  for  the  state  of  an  infinite- dimensional  system 
can  be  computed  only  in  the  simplest  of  situations.  This  is  due  to  the  difficulty  of  solving  a 
partial  differential  equation.  In  general,  it  is  necessary  to  use  a  numerical  approximation  to 
the  state  of  the  system  when  designing  a  controller. 

This  leads  immediately  to  several  questions  regarding  robustness.  Is  the  use  of  finite¬ 
dimensional  approximations  a  valid  technique  when  designing  controllers  for 
infinite-dimensional  systems?  If  so,  which  approximation  methods  can  be  used  to  design 
controllers  which  will  perform  as  designed  when  implemented  on  the  actual  system,  and  how 
high  an  order  is  required? 

A  wide  class  of  projection  methods,  which  includes  mode  truncation  and  finite  element 
Galerkin  techniques,  is  considered  in  this  paper.  It  is  shown  that,  if  the  order  of  an  approxi¬ 
mation  to  a  system  is  sufficiently  high,  then  the  approximating  system  is  stabilizable,  and  a 
controller  designed  using  the  approximation  will  stabilize  the  original  system.  Furthermore, 
the  closed  loop  response  of  the  system  will  be  close  to  the  the  closed  loop  response  of  the 
approximating  system.  This  result  is  independent  of  the  technique  used  for  controller  design. 

Other  researchers  have  studied  convergence  of  solutions  to  the  Riccati  equation  [eg.,  9, 
10,  12],  but  have  assumed  either  uniform  stability  or  stabilizability  of  the  approximating 
systems.  Balas  [1,  2]  assumes  convergence  of  a  stabilizing  feedback  control.  However,  the 
assumptions  on  an  approximation  scheme  required  for  the  results  in  this  paper  to  hold  are 
not  stronger  than  those  typically  imposed  on  any  numerical  technique  used  in  computer 
analysis.  The  key  assumption  is  that  the  original  infinite- dimensional  system  is  stabilizable. 

We  consider  semigroup  control  systems  of  the  following  form: 

x(f)  =  Ax(t)  +  Bu{t), 

y(t)  =  Cx(t),  (1.0) 

x(0)  =  x0  Xo  £  D(A)  C  X. 

A  satisfies  the  Hille-Yosida  theorem  on  the  Hilbert  space  X  and  hence  generates  a  strongly 
continuous  (C^-semigroup  T(t)  on  X.  Control  and  observation  are  assumed  to  be  bounded 
i.e., 

B  e  X)  and  C  E  B(X,  BF). 
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This  formulation  includes  a  wide  class  of  systems;  in  particular  many  problems  concerning 
vibrations  in  structures  with  point  sensing  and  actuation  can  be  placed  within  this  frame¬ 
work. 

The  inner  product  on  X  is  indicated  by  (  ,  )  and  the  corresponding  norm  by,  ||  ||.  Inner 
products  on  different  spaces  are  not  distinguished  as  the  context  indicates  which  space  is 
meant.  We  use  the  following  definition  of  internal  stability. 

Definition  1.1:  The  Co-semigroup  T{t)  is  (7-stable  if  there  exist  constants  M  and  a  >  a 
such  that  ||T(t)||  <  Me~at  for  all  t  >  0. 

The  possibility  of  a  >  0  is  included  to  allow  for  situations  where  the  purpose  of  con¬ 
troller  design  is  for  purposes  other  than  simple  stability  (i.e.,  a  =  0).  For  instance,  while 
vibrating  structures  are  generally  open-loop  stable,  the  settling  time  of  these  systems  is  often 
unacceptably  large. 

Several  definitions  will  be  required. 

Definition  1.2:  The  pair  (A,  B )  is  a-stabilizable  if  there  exists  a  bounded  linear  operator 
K  :  X  — >  Rm  such  that  A  —  BK  generates  a  (7-stable  semigroup. 

Definition  1.3:  The  pair  ( A ,  C)  is  a-detectable  if  there  exists  a  bounded  linear  operator 
F  :  Rp  — >  X  such  that  A  —  FC  generates  a  (7-stable  semigroup. 

Definition  1.4:  The  system  (A,  B,  C)  is  jointly  a-stabilizable /detectable  if  ( A,B )  is 
(7-stabilizable  and  (A,  C)  is  a-detectable. 

Define  the  linear  space  A  to  consist  of  functions  of  the  form 

h  =  ha(t)  +  f;  hx6{t  -  U).  (1.1) 

»=i 

where  ha  €  Z/i[0,  oo)  and  the  sequence  {hi}  is  absolutely  summable.  The  norm  of  a  function 
in  A  is 

11*11  «  f"  \K(t)\dt  +  £;  IM  <  oo.  (1.2) 

J0  i= o 

A  SISO  system  is  said  to  be  Z^-stable  or  bounded-input  bounded-output  (BIBO)  stable 
if  and  only  if  its  impulse  response  belongs  to  the  set  A.  Moreover,  the  gain  of  an  BIBO-stable 
system  is  equal  to  the  norm  of  its  impulse  response  in  the  set  A  [6].  The  symbol  A  is  used 
to  denote  the  set  of  Laplace  transforms  of  distributions  in  A. 

Define  M(  A)  to  be  the  set  of  matrices  with  elements  in  A.  A  MIMO  system  is  Loo-stable 
if  and  only  if  its  impulse  response  belongs  to  Af(A).  The  Loo-gain  of  a  m  x  p  matrix  H  is 
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the  induced  norm  defined  by 


\\H\\a  :=  max  £  II^o  IIa- 
j=i 

Consider  the  common  situation  where  the  system  is  already  stable  and  the  aim  of  con¬ 
troller  design  is  to  improve  its  transient  response.  Let  us  specify  a  real  number  a  >  0,  which 
is  the  minimum  acceptable  stability  margin.  We  define  a  system  to  be  a  —  Loo-stable  if  and 
only  if  its  impulse  response  /(f)  is  of  the  form  (1-1),  and  instead  of  (1.2)  satisfies  the  more 
stringent  conditions 

°°  roo 

V)  \hi\ exp(att)  <  oo,  /  |ft0(i)  exp(crt)|<if  <  oo.  (1-3) 

i=o  Jo 

Let  Ao-  denote  the  set  of  all  distributions  satisfying,  and  let  A,  denote  the  set  of  Laplace 
transforms  of  distributions  in 

We  now  define  precisely  what  is  meant  by  external  stabilization  of  a  system.  Suppose 
P  is  a  given  system,  for  which  we  wish  to  design  a  controller  C ,  arranged  in  the  familiar 
feedback  configuration  shown  in  Figure  1.  It  is  readily  seen  that  the  2x2  transfer  matrix 
H(P,C )  which  maps  the  pair  (ui,u2)  into  the  pair  (ei,  e2)  is  given  by 

'(i  +  pcy1  -Pii  +  cpy1 

H(P,C)  = 

_c(i  +  pc)-'  ( i  +  cp)~1 

The  feedback  system,  or  alternatively  the  pair  (P,  C ),  is  said  to  be  stable  if  each  of  the  four 
elements  in  the  above  matrix  belongs  to  the  set  S  of  stable  transfer  functions.  Definition 
of  S  depends  upon  the  application.  Thus  the  closed  loop  system  is  BIBO-cr-stable  if  and 
only  if  all  four  elements  belong  to  Aa.  It  is  important  to  note  that  the  present  definition  of 
stability  is  symmetric  in  P  and  C.  Thus  P  stabilizes  C  if  and  only  if  C  stabilizes  P. 


u2 


Figure  1.  Feedback  System 
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The  problem  is  to  design  a  controller  for  the  system  (1.0)  which  achieves  such  objectives  as 
external  stability.  Since  this  generally  requires  knowledge  of  the  system  state,  the  difficulty 
in  solving  partial  differential  equation  systems  and  many  other  semigroup  systems  means 
that  an  approximation  to  the  state  must  be  used.  In  the  next  section  we  define  a  class  of 
approximation  techniques  known  variously  as  Galerkin  and  projection  techniques. 

2.  Approximation  Technique 

Suppose  we  have  a  sequence  of  finite-dimensional  subspaces  Xh  with  X ^  C  Xh3  for 
hi  >  h2.  We  assume  that  Xh  C  D(A).  When  finite  elements  are  used  as  a  basis  for  the 
approximating  subspace  Xh,  then  A  should  be  understood  in  its  weak  or  extended  sense  [7, 
20].  Define  PhX  as  the  orthogonal  projection  of  x  £  X  onto  the  finite-dimensional  subspace 
Xh-  The  parameter  h  converges  to  zero  as  the  order  of  the  approximation  is  increased.  Since 
all  norms  on  Xh  are  equivalent,  we  will  assume  that  the  norm  on  Xh  is  that  inherited  from 
X.  Since  Ph  is  an  orthogonal  projection,  it  follows  that 

IM  <  11*11- 

The  approximating  system  is  (Ah,  Bh,  Ch )  where 

Bh  :=  PhB 

Ch  ■=  C\xh 

and  Ah  is  the  bounded  linear  operator  defined  by 

AhXh  =  PhAxh  for  all  xh  £  Xh- 

In  common  with  [3],  the  operators  Ah,  Bh,Ch  and  the  semigroup  Th(t )  generated  by  Ah 
are  only  defined  on  Xh- 

We  make  the  following  assumptions: 

Al)  The  projection  operators  Ph  converge  strongly  to  the  identity  on  the  Hilbert  space 
X.  That  is,  for  all  x  £  X 

lim  || PhX  —  x\\  =  0. 

h  -*0 

A2)  The  approximation  scheme  approximates  the  solution  space  not  the  operator  A,  in 
that  x  £  Xh  implies  Ax  £  Xh  so  that 

AhX  =  Ax  for  x  £  Xh- 

Modal  truncation  and  many  applications  of  the  finite  element  method  satisfy  this  assumption. 
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A3)  A  core  C  of  a  closed  operator  A  is  a  linear  space  contained  in  the  domain  of  A  with 
the  property  that  the  set  of  elements  (x,  Ax),  x  6  C  is  dense  in  the  graph  G(A)  of  the 
operator  A  ([15]  pg,  166).  In  other  words,  if  we  define  an  operator  L  to  be  the  restriction  of 
A  to  C  so  that 

D  (L)  =  C 
Lx  =  Ax,  x  E  C 

then  A  is  the  closure  of  L,  or  equivalently,  G (L)  =  G(A). 

We  will  assume  that  there  exists  a  core  C  for  A  such  that 

lim  \\PhAx  -  AhPhx\\  =  0,  for  all  x  6  C.  (2.1) 

h— *0 

Such  an  approximation  scheme  is  said  to  be  consistent. 

We  will  further  assume  that  the  semigroups  T^t)  generated  by  A/,  are  uniformly  bounded, 
that  is,  there  exist  real  numbers  H,  M  >  1  and  u  such  that 

||Tfc(t)||  <  Meut  for  all  h<H.  (2.2) 

(Uniform  boundedness  of  the  approximate  semigroups  is  generally  referred  to  as  “stability” 
in  the  numerical  analysis  literature.)  Consistency  and  uniform  boundedness  are  required  for 
convergence  of  the  approximation  [e.g.,  16,  18,  20]  i.e.,  for  all  e  >  0,  i  >  0  and  for  all  x  6  X 
there  exists  H  such  that 

|| PhT(r)x  —  Th(r)Phx\\  <  e  for  all  r  £  [0,1]  and  h  <  H. 

Assumptions  (A1)-(A3)  are  satisfied  by  typical  approximation  methods  such  as  finite 
element  techniques  and  modal  truncation. 

In  subsequent  sections  the  original  semigroup  control  system  will  be  assumed  stabilizable 
and/or  detectable.  However,  no  assumptions  about  the  stabilizability  or  detectability  of 
the  approximations  are  made.  In  the  next  section  it  is  shown  that  stabilizability  of  the 
approximating  systems  can  be  deduced  from  stabilizability  of  the  original  semigroup  control 
system. 

3.  Stabilizability  of  Approximating  Systems 

We  require  the  following  theorem  which  guarantees  the  existence  of  a  solution  to  the 
infinite-dimensional  Riccati  equation  if  the  pair  (A,B)  is  stabilizable. 

Theorem  3.1:  [17]  Let  A  and  B  be  as  defined  above  and  let  L,Q  be  self-adjoint  and 
positive  definite  bounded  linear  operators.  If  there  exists  a  solution  u*  to  the  optimal  control 
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problem  with  the  quadratic  cost  function 

J(u)  :=  f  [(Lx(t),  x(t))  +  (Qu(t),u(t))}dt  (3.1) 

Jo 

which  satisfies  the  condition  that  J(u*)  be  finite,  then  there  exists  an  unique  self-adjoint  and 
non-negative  solution  IT  £  L(Aj  A)  to  the  infinite-dimensional  Riccati  equation 

A*U  +  UA-UBQ-1B*n  +  L  =  0  (3.2) 

where  0  is  the  zero  operator.  Let  K  =  Q_1B*II;  then  A  —  BK  generates  an  exponentially 
stable  Co-semigroup.  □ 

Clearly  an  admissible  solution  exists  if  and  only  if  (A,  B )  is  stabilizable. 

Theorem  3.2:  Assume  we  have  a  semigroup  control  system  (2.0)  and  an  approximation 
scheme  satisfying  assumptions  (A1)-(A3)  in  Section  2.  Assume  that  the  system  is  stabiliz¬ 
able.  Then,  for  sufficiently  small  h,  the  approximating  systems  are  stabilizable  i.e.,  there 
exists  an  operator  Kh  such  that  Ah  —  BhKh  generates  an  exponentially  stable  semigroup  on 

xh. 

Proof:  For  each  approximating  subspace,  define  a  self-adjoint  operator  Uh  ■  Xh  —>  Xh 
as 

Bh  :=  PhUPh 

where  II  is  the  solution  to  the  Riccati  equation  (3.2).  It  will  be  shown  that,  for  small  enough 
h,  II/i  solves  the  finite-dimensional  Riccati  equation 

A^IIfc  +  rU-A/,  -  II hBhQ~l B'hRh  +  Lh  —  0^ 

where  0^  is  the  zero  operator  in  B(A),,  Xh)  and  Lh  is  a  self-adjoint  positive  definite  operator 
(as  yet  undetermined)  in  B(A'/,,  Xh). 

Define  the  operator  R{Lh)  :  Xh  — >  Xh  by 

r(lh)  =  A*hnh  +  nhAh  -  n +  lh 

where  Lh  £  B(A\,  Xh)  is  as  yet  unspecified.  Let  us  calculate  a  representation  lor  R(Lh) 
with  respect  to  some  basis  {<£,}"=1  for  Xh-  Using  (3.2),  we  have 

R{Lh)h)  =  (*,  QWh)  +  (*,  Mi)  -  (*,  L<f>j) 

where 

Q(h)  :=  UBQ-'B'Il  -  PhUPhBQ-'B'PhllPh- 
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Now  define  Lh  ■  Xh  —*  X h  as 


Lhxh  =  Lxh  -  Q(h)xh  for  all  xh  £  Xh. 

Then  R(Lh)  =  0  on  H(Xh,Xh)  since 

( <t>i ,  R(Lh)<t>j )  =  0  for  all  <f>t,  <f>r 

In  order  to  show  that  IT*,  solves  a  Riccati  equation  for  (Ah,  Bh )  for  small  enough  A, 
it  remains  only  to  prove  that  Lh  is  self-adjoint  and  positive  definite  for  small  h.  Now, 
lim/,_o  || (I  —  -P/»)-0||  —  0  since  B  operates  on  a  finite-dimensional  space.  It  follows  that  [22, 
Theorem  4.14] 

lim||g(A)||  =  0. 

The  fact  that  Lh  is  self  adjoint  follows  easily  from  the  definition  of  Q(h).  Since  L  is  positive 
definite  and  self-adjoint  we  can  choose  H  so  that  for  all  h  <  H 

(x/i,  Lh%h>)  ^  C  (xh,  Xh) 

for  some  c  >  0.  Thus  Lh  is  self-adjoint  and  positive  definite  for  all  h  <  H,  and  II/,  PhLLPh 
solves  the  Riccati  equation  for  ( Ah,Bh )  with  weights  Q  and  Lh  if  h  is  sufficiently  small. 
Defining  Kh  Q~lB^Hh,  Ah  —  BhKh  generates  a  stable  semigroup  on  Xh  for  h  <  H.  □ 

Banks  and  Kunisch  [3]  derived  a  similar  result  for  parabolic  systems  using  different  as¬ 
sumptions.  Rosen  [19]  considered  the  problem  of  convergence  of  solutions  to  approximating 
Riccati  equations.  However  he  considers  a  much  smaller  class  of  problems:  his  proof  applies 
to  cases  where  the  open  loop  system  is  stable  and  II  is  a  Hilbert-Schmidt  operator.  Although 
the  result  given  here  is  weaker  in  that  convergence  of  11^  to  II  in  an  operator  norm  is  not 
shown,  the  theorem  is  considerably  more  general. 

The  following  corollary  is  a  consequence  of  the  two  results  above. 

Corollary  3.2a:  Assume  that  ( A,B )  is  cr-stabilizable  and  let  n  be  the  solution  to 
the  Riccati  equation  for  (A  +  al,  B)  with  weights  Q  and  R  as  in  equation  (3.2).  Define 
Ilh  :=  PhUPh  and  let  Kh  =  Then  for  sufficiently  small  h,  Ah  —  BhKh  generates  a 

cr-exponentially  stable  semigroup  on  Xh-  □ 

The  results  above  for  cr-stabilizable  systems  have  obvious  extensions  to  a-detectable 
systems  since  ( A,C )  is  a-detectable  if  and  only  if  (A*,C*)  is  a-stabilizable. 

Theorem  3.3:  Let  (A,  B,  C )  be  a  semigroup  control  system  and  suppose  we  have  an  ap¬ 
proximation  scheme  which  satisfies  assumptions  (A1)-(A3)  of  Section  2.  Assume  furthermore 
that  ( A,B )  is  a-stabilizable.  Choose  II  £  B(X,X)  such  that  II  solves  the  Riccati  equation 
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(3.2).  Define  K  Q~1B*Tl  and  let  S(t )  to  be  the  semigroup  generated  by  A  —  BK .  Simi¬ 
larly,  define  Kh  ■=  Q_1  B^PhJlPh  and  let  Sh(t )  be  the  semigroup  generated  by  Ah  —  BhKh- 
Then  S/,(f )  converges  to  S(t)  uniformly  on  bounded  intervals;  i.e.,  for  all  e  >  0,t  >  0  and 
for  all  x  £  X  there  exists  H  such  that 

|| PhS(r)x  —  5^(r)P^x||  <  t  for  all  r  £  [0,  t]  and  h  <  H. 

Proof: 

\\BhKh\\  <  \\BQ-'B'\\  ||n||  :=  c. 

Using  (2.2)  and  [18,  Theorem  3.1.1]  it  follows  that  there  exist  real  numbers  H,M  >  1  and 
o u  such  that 

\\Sh{t)\\  <  Me(“+Wcl‘  for  all  t  >  0,  for  all  h  <  H.  (3.3) 

Define  Aa  A  —  BK  and  Aho  :=  Ah  —  BhKh-  Let  C  be  the  core  for  A  such  that  the 
consistency  condition  (2.1)  is  satisfied.  Since  BK  is  a  bounded  operator,  C  is  also  a  core  for 
A0.  The  resolvent  of  (XI  —  Aa )  is  bounded  for  Re(A)  >  - a ,  and  so  (XI  -  j40)C  is  dense  in 
X  for  ReA  >  -  a.  It  follows  that  for  every  x  £  C,  we  have 

Jim  \\AhoPhX  -  PhAoX ||  <  Jim  \\PhBQ~1  B*(PhUPh  -  II)x||  =  0. 

h — »0  h—*0 

In  order  to  show  convergence  of  {S^}  to  S  it  necessary  to  introduce  a  semigroup  which  is 
defined  on  all  of  X.  Define  Aho  —  AhoPh-  Since  Aho  is  a  bounded  linear  operator  it  generates 
the  Co-semigroup  Sh(t)  defined  by 

=  E  T-^l  =  £  =  ShWPh. 

k=0  K-  k= 0  K- 

Since  ||5/l(t)||  <  ||5h(t)||,  Sh(t)  also  satisfies  (3.3).  Now,  for  x  £  C,  we  have 

lim  ||  AhoX  -  Ax ||  <  lim  ||  AhoPhx  -  PhA0x\\  +  ||Pfci4„x  -  i4„x||. 

The  first  term  was  shown  above  to  tend  to  zero  .  The  second  term  converges  to  zero  since 
Ph  converges  strongly  to  the  identity.  It  follows  that 

lim  sup  ||S)i(t)x  —  S'(t)x||  =  0  for  all  x  £  X,  for  all  r  >  0 

0  0<t<T 

by  the  Trotter-Kato  Approximation  Theorem  [18,  Theorem  4.5].  Therefore,  for  all  x  £ 

X,r  >  0, 

lim  sup  ||5^(t)Phx  -  PhS(t)x ||  =  0. 

h—0  0<t<T 


The  following  theorem  is  now  immediate. 

Theorem  3.4:  Suppose  that  the  system  ( A,B,C )  is  stabilizable  and  that  assumptions 
(A1)-(A3)  hold.  Assume  also  that  the  approximation  scheme  also  satisfies  the  additional 
condition  (A4): 

lim  sup  || Tfr(t)PhX  —  PhT*(t)x\\  —  0. 

h~>0  o  <t<r 

Then  for  sufficiently  small  h,  the  Riccati  equation 

A*hrh  +  7 rhAh  -  7T hBhQ~xB*h-Kh  +  PhLPh  =  Qh  (3.4) 

has  a  unique  non-negative  self-adjoint  solution  7r^  .  The  sequence  of  operators  7 converge 
strongly  to  the  solution  to  the  exact  equation  (3.2).  Define  the  optimal  feedback  operators 
K) j  and  K  :=  R~1B II;  the  semigroup  generated  by  A /,  —  B^Kh.  converges  to 

the  semigroup  generated  by  A  —  BK. 

Proof:  By  assumption,  hypotheses  (H1)-(H2)  in  [12]  are  satisfied.  Since  Theorem  3.3 
implies  uniform  stabilizability/detectability  it  follows  that  the  sequence  of  solutions  to  the 
approximate  Riccati  equations  are  bounded  by  some  real  Mj  [12].  Since  {7^}  £  B(H,H) 
the  sequence  contains  a  weakly  convergent  subsequence.  It  is  clear  that  the  limit  of  this 
subsequence  is  also  self-adjoint  and  non-negative.  In  fact,  the  limit  is  the  solution  to  the 
exact  Riccati  equation,  II  [Theorem  6.7,  8].  A  proof  identical  to  that  of  [Theorem  6.9,  8] 
yields  the  remainder  of  the  result.  □ 

The  fourth  assumption  can  easily  be  shown  to  hold  for  a  wide  class  of  problems.  How¬ 
ever,  while  assumptions  (A1)-(A3)  are  almost  always  satisfied  by  a  numerical  scheme  which 
is  satisfactory  for  approximation  purposes,  (A4)  may  not  be  satisfied.  This  can  lead  to 
unsatisfactory  controller  design.  An  example  of  this  situation  is  discussed  in  [4], 

4.  Convergence  in  the  Graph  Topology 

The  graph  and  gap  topologies  [22,  23]  arose  from  a  need  to  define  convergence  of  possibly 
unbounded  operators.  The  basic  idea  was  outlined  in  Kato  [15].  If  M,N  are  closed  linear 
subspaces  of  a  Banach  space  X  then  S(M,N )  is  the  smallest  number  6  such  that 

6  >  dist(u,  N)  =  in f  ||u  —  v||  for  all  u  £  M.  (4-1) 

The  gap  between  M  and  N  is  defined  by 

6{M,N )  :=  ma x{6(M,  N),6{N,M)). 
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If  S  and  T  are  closed  operators  from  X  to  Y  then  the  gap  between  their  graphs  as  closed 
subspaces  of  X  x  Y  is  well  defined.  The  difficulty  in  defining  a  gap  for  unstable  systems  is 
that  the  plant  may  not  be  represented  by  a  mapping  between  two  Banach  spaces.  However, 
some  bounded  inputs  u,  will  be  mapped  to  bounded  outputs  y.  Let  U  be  the  set  of  stable 
inputs  and  outputs  and  let  P  be  the  impulse  response  of  a  scalar  plant  .  We  will  assume 
that  U  is  a  Banach  space.  The  symbol  denotes  convolution.  Define 

Du  :=  {u  £  U  such  that  y  £  U} 


and  the  operator 


P\j  *  u  :=  P  *u, 


with  domain  Du-  Then  the  graph  of  i\j  is 


G(Pv)  =  {(it,!/),  y  =  P  *u,u  £  Du}. 

It  is  easy  to  show  that  G(P\j)  is  closed  if  P  is  a  closed  operator.  The  graph  of  a  plant  is 
now  defined  to  be  the  set  of  bounded  (in  the  appropriate  sense)  input-output  pairs,  be., 

Gu (P)  =  G(Pv)  =  {(«,»)  £  U"“+',  y  =  P  *  u}. 

In  the  case  of  BIBO  stability, 

G„(/>)  =  {(u,so£irp,}. 

Define  the  gap  between  operators  as 

d(P1,P2)  =  ^(Gu(Fi),Gu(P2)). 

With  this  metric  the  space  of  closed  operators  from  U  to  U,  C(U,U)  becomes  a  metric 
space.  Convergence  of  a  sequence  Pn  — ►  P  is  defined  by  d(Pn,  P)  — ►  0.  If  a  sequence  of 
plants  {P,}  are  stable,  then  convergence  in  norm  to  P  is  equivalent  to  convergence  in  gap. 
The  topology  induced  by  the  metric  d( Pi ,  P2 )  is  thus  an  appropriate  generalization  of  the 
uniform  topology  to  unstable  plants. 

Note  that  the  above  definition  of  “stable  graphs”  includes  the  graphs  of  non-causal  plants. 
A  more  useful  topology  for  our  purposes,  which  considers  only  linear  causal  plants,  can  be 
obtained  through  the  frequency  domain  approach  of  Vidyasagar  [22].  Let  the  set  of  stable 
plants  S  be  A,  P  a  plant  with  a  right  coprime  factorization  (N,  D)  over  A.  Define  G^P) 
to  be  the  Laplace  transform  of  elements  of  Goc(P).  Then 

G00(P):={(u,AP-1u),u£  LZ} 
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or 


G  s(P)^{(Dz,Nz),zeLZ}. 


The  graph  topology  for  unstable  plants  is  defined  in  [22]  using  right  coprime  factorizations 
(r.c.f.’s).  Essentially,  a  basic  neighborhood  in  the  graph  topo’ogy  of  a  plant  with  a  r.c.f. 
(N0,  Da)  is  the  set  of  all  plants  with  a  r.c.f.  (N,D)  where  N  is  close  to  N0  in  norm  and 
D  is  close  to  Da.  Define  S  to  be  the  usual  algebra  of  stable  rational  functions.  For  finite¬ 
dimensional  plants  the  graph  topology  is  the  topology  induced  on  the  quotient  space  S/S 
by  the  set  of  stable  functions,  S.  Further  details  can  be  found  in  [22].  Zhu  [24]  presents  a 
generalization  of  the  graph  and  gap  topologies  to  subsets  of  general  quotient  fields.  Here, 
we  consider  only  a  particular  set. 

Callier  and  Desoer  [5]  define  a  quotient  algebra  of  transfer  functions  of  distributed  systems 
which  is  useful  for  studying  control  of  semigroup  systems. 

We  say  that  a  transfer  function  /  in  A  is  bounded  away  from  zero  at  infinity  if  for  every 
sequence  {s*}  with  Re(s)  >  0  and 

1**1  oo 


we  have 


Define  the  sets 


lim  inf  |A(sfc)|  >  0. 

k~*  oo 

A_  =  {h  £  Aa  for  some  a  >  0} 

A_i0O  =  {h  £  A _,/  bounded  away  from  zero  at  infinity} 


and  the  corresponding  quotient  algebra 


B  =  {h  =  a/b,a  £  A_,6  £  A_)CO}. 

The  assumption  that  b  €  A_iOC  implies  that  functions  in  B  possess  only  a  finite  number  of 
unstable  poles.  The  extension  to  the  multivariable  case  is  straightforward,  and  we  define 
M(B)  to  be  the  set  of  transfer  function  matrices  with  entries  in  B. 

An  extension  which  restricts  the  region  of  stability  to  Re(s)  <  —  a  is  also  straightforward. 
The  algebra  A  is  replaced  by  Ka  in  the  above  definitions  and 

Ba  =  (A  =  a/b,  a  £  Aa>_,  b  £  Ar.oo}- 


Jacobson  [13]  and  Jacobson  and  Nett  [14]  have  shown  that  semigroup  control  systems 
which  are  o-stabilizable  and/or  cr-detectable  have  transfer  functions  in  M(Bff).  This  frame¬ 
work  is  thus  appropriate  for  studying  semigroup  control  systems  using  frequency  domain 
methods.  Furthermore,  every  function  P  in  M{ Bff)  has  both  a  left-  and  a  right-coprime 
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factorization  over  M( Aa).  Thus  Vidyasagar’s  graph  topology  is  well-defined  for  semigroup 
control  systems  providing  they  are  a-stabilizable  and/or  -detectable. 

Theorem  4.1:  Suppose  Fh  is  a  sequence  in  M(Ba),  and  that  F  G  M(B<7).  Then  the 
following  statements  are  equivalent. 

(i)  {F*}  converges  to  F  in  the  graph  topology. 

(ii)  There  exist  a  r.c.f  ( N,D )  of  F ,  and  a  sequence  of  r.c.f.’s  (Nh,Dh)  of  Fh  such  that 
Nh  — ►  N  and  D h  — ►  D  in  M(Aa) 

(iii)  There  exist  a  l.c.f.  ( N,D )  of  F,  and  a  sequence  of  l.c.f.’s  (Nh,Dh)  of  Fh  such  that 
Nh  — >  N  and  Dh  — ♦  D  in  M(AtT)  . 

Proof:  Since  all  transfer  functions  in  M( B^.)  have  left  and  right  coprime  factorizations 
[5],  Vidyasagar’s  proof  [21]  for  rational  functions  that  convergence  in  the  graph  topology  is 
equivalent  to  convergence  of  coprime  factorizations  extends  in  a  straightforward  manner  to 
transfer  functions  in  M(Btr).  □ 

The  following  result  can  be  loosely  paraphrased  as  follows:  A  family  of  plants  Fh  can  be 
robustly  stabilized  by  a  compensator  C  which  stabilizes  some  nominal  plant  F  if  and  only 
if  Fh  converges  to  F  in  the  graph  topology.  Furthermore,  the  closed  loop  response  of  the 
feedback  pair  ( Fh,C )  converges  to  that  of  ( F,C ).  In  other  words,  the  graph  topology  is  the 
weakest  topology  in  which  feedback  stability  is  robust. 

Theorem  4.2:  Let  Fh  be  a  sequence  of  plants  in  . 

(i)  Suppose  Fh  converges  tc  F  €  Af( Ba)  in  the  graph  topology.  Let  C  G  M( B^.)  stabilize 
F.  Then  there  exists  an  H  such  that  C  stabilizes  Fh,  for  all  h  <  H,  and  moreover,  the 
closed  loop  transfer  matrix  H(Fh,C)  converges  to  H(F,C)  in  M(A<r). 

(ii)  Conversely,  suppose  there  exists  a  C  G  M(Ba)  which  stabilizes  Fh  for  all  h  <  H,  and 
that  H(Fh,  C)  converges  to  H(F,  C).  Then  {F/,}  converges  to  F  in  the  graph  topology. 

Proof:  Identical  to  that  in  [21].  □ 

This  has  obvious  implications  for  controller  design  using  approximations.  Failure  of  a 
sequence  of  approximations  to  converge  in  the  graph  topology  implies  that  at  least  one  of 
the  following  must  occur  (1)  simultaneous  stabilization  of  Fh  and  F  is  not  possible  for  all  h 
sufficiently  small  or  (2)  convergence  of  the  closed  loop  response  H(Fh,  C )  does  not  occur. 
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Theorem  4.3:  Let  (A,  5,(7)  be  a  o-stabilizable/detectable  semigroup  control  system, 
and  suppose  we  have  a  sequence  of  approximations  (.A/,,  Bh,  Ch)  satisfying  assumptions  (Al)- 
(A3)  of  Section  2.  If  the  system  is  cr-stabilizable  and  C  is  any  controller  which  stabilizes  P, 
then  for  sufficiently  small  H,  (P/,,  C )  is  cr-stable  for  all  h  <  H.  Furthermore,  the  closed  loop 
response  of  (P/,,  C )  converges  in  norm  to  that  of  (P,  C). 

Proof:  The  theorem  will  be  proven  by  using  Theorem  3.3  to  show  that  the  approximate 
systems  (A/,,  Bh,  Ch)  converge  to  (A,  B ,  C)  in  the  graph  topology,  and  then  applying  Theorem 
4.2. 

Let  feedback  operators  K  and  Kh  be  as  defined  above  so  that  A^  :=  A/,  — 5/,  A),  generates 
a  cr-stable  semigroup  Sh(t)  for  sufficiently  small  h  and  A0  :=  A  —  BK  generates  a  cr-stable 
semigroup  S(t).  For  small  h,  right  cr-coprime  factorizations  for  the  approximate  systems  are 
given  by 

M3)  —  '■  Aho)Bh , 

Dk(- s)  =  I  —  KhR{s  :  Aho)Bh 

where  P(s  :  Aho)  indicates  the  resolvent  of  A^>  [21].  A  right  coprime  factorization  for  the 
original  system  is  [13] 

N(s)  =  CR{s  :  Aa)B, 

D(s)  =  /  -  KR(s  :  A0)B. 

The  result  follows  if  it  can  be  shown  that  Nh{t)  —>  N(t)  and  Dh(t)  — >  D(t)  in  the  norm 
on  A*.  Convergence  of  the  numerators  is  proven  first. 

Since  {5^(0}  converge  uniformly  on  bounded  intervals  to  S(t),  there  exist  H  and  M  >  1 
such  that  for  some  7  >  a 

||5(0II  ^  Me-7t  for  all  t  >  0, 

and 

||5*(0||  <  Me’7'  for  aU  t  >  0,  for  all  h  <  H.  (4.2) 

Therefore,  for  any  e  >  0  there  exists  r  such  that 

jr  exp((r0||Af(0  -  Nh(t)\\dt  <  |.  (4.3) 

Also, 

sup  \\N(t)  -  JVK(i)||  <  sup  IICll  II S(t)B  -  Sh(()«II-  (->.4) 

0<«<T  0<t<T 

It  follows  from  Theorem  3.3  that  this  also  approaches  zero.  Combining  statements  (4.3)  and 
(4 A)  we  see  that 

lim  f  exp(ai)|]iV(0  -  Ah(t)P  =  0 

h<0  J 0 
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and  so  N ^  converges  to  N  on  A,. 

The  proof  showing  convergence  of  the  denominators  is  somewhat  lengthier  and  requires 
use  of  a  lemma  in  [16].  Defining 

Fh(t)  :=  (Ph n  -  n )Sh(t)PhB, 


we  have 

Dh(t)  -  D(t)  =  Q~lB-Fh(t)  +  Q^BmiS^PnB  -  S(t)B}.  (4.5) 

lim  ||  [°°  e~XtFh(t)dt\\  <  lim  ||(PJI  -  II)  /°°  e~XtSh(i)Bdt\\ 
h~*  0  Jo  h-rO  Jo 

<  Wm\\{PhU-U)R(\-,Aho)B\\ 

h — >0 

<  lim  |!PhII||(jR(A;  Aho)B  —  R(\;  A0)B\\ 

h—*0 

+  ||(P*-/)n*(A;A0)B|| 

+  PH  \\(R(X-,A0)B-R{X-Aho)B\\ 

=  0, 

since  Ph.  converges  strongly  to  the  identity  and 

lim  || J2(A;  A^B  -  R( A;  A0)B\\  =  0 

h — *0 

for  all  A  [18,  Theorem  4.2].  The  sequence  will  now  be  shown  to  be  equicontinuous  . 

The  core  C  is  dense  in  D(A)  which  is  dense  in  X.  Let  {ej}"=1  be  an  orthonormal  basis 
for  Rm  and  for  any  6  >  0  define  a  bounded  operator  Gf  :  RF  — >  X  by 

Gg(ei)  =  X{  where  x^  G  C  and  UBe^  —  xj|  <  S/m. 

Then  G$  is  a  bounded  linear  operator  with  range  contained  in  C  and  with 

\\GS  -  B\\  <  6. 

Now  is  bounded  by  some  real  Mi  since  it  converges  to  ||.A0G$||.  Therefore, 

||f»(t  +  s)  -  f\(f) II  <  || /yi  -  nil  II  !*  Sh(t  +  ^A^BduW 

<  liny  ||jS^(i  +  u)ifeG,ju||+ 

<  ||n||  |||5h(t  +  S)  -  S„(t)](B  -  G,)||. 
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Let  M  be  the  constant  bounding  Sh(t)  for  small  h  (Equation  5.2).  Then 


||*fc(*  +  6)-  Fh(t) ||  <  \\U\\[SMMl  +  2 M6\ 
and  for  sufficiently  small  H  and  any  e  >  0  we  can  choose  6  so  that 

|| Fh(t  +  6)-  Fh(OII  <  e  for  all  h  <  H. 

Thus  is  an  equicontinuous  family  and  since  i'k(i)  is  also  clearly  uniformly  bounded  for 
all  t  and  h,  it  follows  from  [16,  Lemma  2.11]  that 

lim  sup  || Fh(t)  -  F(0H  =  0. 

0  <t<r 

Since  II  and  Q~lB *  are  bounded  operators  it  follows  from  Equation  (4.5)  and  Theorem  3.3 
that 

lim  sup  || Dh(t)  -  D(t) ||  =  0. 

0  0<t<T 

Finally,  the  remainder  of  the  proof  of  convergence  of  £)/,(£ )  to  D(t)  follows  as  for  the  numer¬ 
ators  (4.3). 

Convergence  in  the  graph  topology  follows  from  Theorem  4.1.  The  conclusions  of  the 
theorem  then  follow  from  Theorem  4.2.  □ 

The  result  is  independent  of  the  technique  used  for  stabilization,  and  does  not  require 
any  assumptions  on  the  stabilizability  or  detectability  of  the  approximations.  For  instance, 
let  ( Nh,Dh )  be  a  cr-stable  r.c.f.  of  an  approximating  plant  and  choose  X,  Y  G  M( Aa)  such 
that 

XNh  +  YDh  =  /. 

Then  if  the  original  plant  is  close  enough  to  the  approximation  so  that,  for  some  r.c.f  of  the 
original  plant,  ( N,D ), 

N-Nh 
D-Dh 

then,  any  controller  designed  using  the  approximating  plant  will  stabilize  the  original  plant. 
Furthermore,  the  closed-loop  responses  will  be  close  in  norm,  and  the  system  will  be  internally 
er-stable.  The  difficult  part,  of  course,  is  determining  whether  (4.6)  is  satisfied. 

5.  Conclusions 

The  above  results  validate  the  use  of  Galerkin  type  approximations  in  controller  design 
for  infinite-dimensional  systems,  regardless  of  the  control  technique ,  and  with  very  weak 
assumptions  on  the  approximation  method.  The  key  assumption  is  that  of  stabilizability/ 
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detectability  of  the  semigroup  control  system  which  supports  the  work  of  Jacobson  and  Nett 
[14]  indicating  the  importance  of  this  property  in  control  of  infinite- dimensional  systems. 

A  high  order  Galerkin  model  can  be  chosen  and  used  to  design  a  controller  using  lumped 
parameter  techniques.  If  the  objective  is  optimal  state  feedback,  then  convergence  of  the 
Kalman  gains  can  be  used  as  a  criterion.  Otherwise,  robustness  of  the  approximate  plant  with 
respect  to  the  original  infinite-dimensional  system  can  be  estimated  before  any  controller  is 
implemented  by  either  (1)  comparison  of  coprime  factors  (if  available)  or  (2)  determination 
that  the  sequence  of  approximations  has  converged  satisfactorily.  Several  possible  criteria  are 
given  in  [21]  and  [11].  An  open  problem  is  determining  when  the  approximation  order  is  high 
enough  for  satisfactory  controller  design.  Since  exact  transfer  functions  will  be  unavailable 
for  most  practical  problems,  a  useful  result  would  be  the  rate  of  convergence  of  the  coprime 
factors,  compared  to  that  of  the  numerical  scheme. 

A  subsequent  paper  will  extend  the  results  in  this  paper  to  discrete  time  approximation 
schemes. 
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